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CN ■ Abstract 

, We develop a new variant of the wave-packet analysis and solve the tunneling 

^ ' time problem for one particle. Our approach suggests an individual asymptotic de- 

scription of the quantum subensembles of transmitted and reflected particles both at 
the final and initial stage of tunneling. We find the initial states of both subensem- 
pLiI bles, which are non-orthogonal. The latter reflects ultimately the fact that at the 

initial stage of tunneling it is impossible to predict whether a particle will be trans- 
mitted through or reflected off the barrier. At the same time, in this case, one 
^ I can say about the to-be-transmitted and to-be-reflected subensembles of particles. 

We show that before the interaction of the incident packet with the barrier both 
the number of particles and the expectation value of the particle's momentum are 
^ , constant for each subensemble. Besides, these asymptotic quantities coincide with 

^ those of the corresponding scattered subensemble. On the basis of this formalism 

we define individual delay times for both scattering channels. Taking into account 
the spreading of wave packets, we define the low bound of the scattering time to 
describe the whole quantum ensemble of particles. All three characteristic times 
are derived in terms of the expectation values of the position and momentum op- 
erators. The condition which must be fulfilled for completed scattering events is 
derived. We propose also the scheme of a gedanken experiment that allows one to 
verify our approach. 
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Introduction 



The tunneling time problem is perhaps one of the most long-lived and involved problems 
of quantum mechanics. Although there exists a variety of proposals (see, for example, 
reviews [1, 2, 3]) to define characteristic times for a tunneling particle, this question have 
been continued to be controversial. The first attempts (see, for example, [1, 4, 5, 6, 7]) 
to solve this problem, on the basis of the one-dimcnsional Schrodingcr equation (OSE), 
were reduced to the naive analysis of the relative motion of the transmitted (or reflected) 
and incident wave packets. Such an approach is usually referred to as the standard wave- 
packet analysis (SWPA). The main result of the SWPA is the development of concepts 
of the phase and delay times to describe particles with a definite momentum. However, 
even in this particular case these times proved to be ill-defined. Their rigorous justifi- 
cation leads to the difficulties which remain to be overcome (see [1, 2]). But the most 
serious difficulties in the SWPA arose in attempt to define tunneling times for wave pack- 
ets whose width is comparable with the width of a potential barrier. Many authors (see, 
for instance, [1, 2]) pointed out that some properties of tunneling such wave packets (the 
dependence of the phase and delay times on the initial distance between the incident wave 
packet and the potential barrier; the acceleration of the transmitted wave packet by an 
opaque potential barrier) seem to be paradoxical. Due to this fact the SWPA has been 
considered (see [2, 3]) as the most unpromising approach to the solution of the tunneling 
time problem. As a result, all the significant papers published on this problem for the last 
ten years have been fulfilled out of the framework of the SWPA (see reviews [2, 3] and 
also [8, 9]). However, as is seen from [2, 3], the alternative approaches themselves meet 
with the serious difficulties. The basic assumptions (to prescribe the Feynman, Bohmian 
or Wigner trajectories to a particle; to incorporate in quantum mechanics the conditional 
probabilities which may be negative by value; to introduce the special "quantum" clocks 
for timing the particle's motion; to admit the averaging procedures (for example, averag- 
ing over time intervals) to differ from the averaging over a quantum ensemble; to admit 
making use, in addition to parameter t, of operators with the dimension of time) under- 
lying these approaches are either purely speculative or questionable in many respects. In 
any case, such assumptions lead the alternative approaches out of the framework of the 
conventional quantum mechanics and, as a consequence, they themselves need an inde- 
pendent justification. A detailed analysis of this problem can be found in [2, 3, 10] (see 
also [11]). 

We have also to add that none of them explains the mentioned above paradoxical 
behavior of wave packets. The point is that such an explanation have not merely been a 
purpose of the alternative approaches. As was said by Landauer and Martin [2], this is 
"...a responsibility of wave-packet analysis proponents...". However, in our opinion, such 
a viewpoint is at least strange because the description of the particle tunneling in terms of 
wave packets is inherent to quantum theory. So that an exhuastive interpretation of the 
wave-packet dynamics in tunneling is a task not only for the proponents of the SWPA. All 
approaches whose aim is to solve the tunneling time problem must explain the behavior 
of wave packets in tunneling. 

Taking into account the above, we see a necessity to return to the SWPA, and clarify 
the reasons having led to the failure of the earlier approaches in solving the given problem. 
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We begin our analysis with the remark that there is really the case in the SWPA when 
no problems arise in timing the particle's motion. We have in mind the motion of an 
everywhere free particle. Namely the timing procedure for a free particle was used in 
defining the phase times (see [1, 4, 5, 6]). The possibility to use such a procedure in the 
tunneling problem is based on the fact that a tunneling particle moves freely before and 
after its interaction with the potential barrier. 

To avoid some mistakes in the more complex, scattering case, one has to dwell shortly 
on the free-particle timing procedure. A simple analysis of classical dynamics of particles 
shows that it is important to distinguish two kinds of characteristic times: 1) the time of 
arrival of a free particle at some given point, and 2) the duration of passing the particle 
through some given spatial interval. Defining the first quantity is based in quantum 
mechanics on the fact that the expectation values of the velocity and position of a free 
particle coincide with the corresponding characteristics of the "center of mass" (CM) of a 
wave packet to describe the particle state. As is known, in this case the motion of the CM 
is described by Newton equations for a free classical particle. Thus, in order to define the 
(mean) time of arrival of a quantum particle at a given point, one needs simply to follow 
the CM of the corresponding wave packet. This concept will be used here in defining 
delay times. 

A more distinct difference between the classical and quantum dynamics of free particles 
arises in determining the characteristic times of the second kind. So, the time a particle 
spends within some spatial region as well as the time the particle spends to traverse this 
region are the same classically. However, this is not the case in quantum mechanics. For 
example, a quantum particle spends no time in the null spatial interval. But there is 
a non-zero temporal interval when the probability for a particle to traverse this region 
is large enough. To define the duration of the physical process, one has to follow at 
least the fronts of a wave packet, rather than its CM. For this purpose, in addition to 
the expectation value of the position operator at some instant t, one needs to know its 
mean-square deviation. This concept will be used below in defining the low bound of the 
scattering time. 

The free-motion case reveals the following two aspects of timing a quantum dynamics 
of particles. Firstly, timing the motion of a quantum particle can be thought only as 
timing the evolution of different (initial or central) moments of the position (or another 
Hermitian) operator. This is related ultimately to the fact that quantum mechanics can 
predict only expectation values of physical quantities. Secondly, there is no necessity 
to construct special "quantum clocks" for timing the motion of a free quantum particle. 
Characteristic times for such a particle (in other words, for the CM or fronts of a free 
wave packet) can be measured with the help of the same clocks that used in keeping track 
of a classical free particle. 

It is obvious that the above peculiarities of timing a free particle must be true also 
for a scattering particle when it moves beyond the scattering region (this suggests that 
the scattering event is a completed scattering). However, in this case a new aspect of 
timing arises: while a particle interacts with the barrier, to time its motion is impossible. 
There are two main reasons for this. Firstly, while a wave packet interacts with the 
barrier, the position and velocity of its CM say nothing about the expectation values of 
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the position and velocity of a particle: because of the interference, in this case they can 
be unambiguously related neither to an incident nor transmitted, nor reflected particle. 
Secondly, any measurement in the barrier region contradicts the conventional setting of 
the tunneling problem. 

As is seen, the free-particle timing procedure itself is simple enough. However, the 
application of this procedure to a scattering particle, in the case when the wave packet 
width is comparable with the barrier width, leads to the ill-defined phase times (see, for 
example, [1, 2, 4, 12, 13]. And the formalism based on this procedure does not provide a 
proper explanation to the behavior of wave packets in tunneling. As is known, a numerical 
modelling of the tunneling process for wide (in x-space) wave packets shows that, though 
the interaction of a particle with the barrier is elastic, the CMs of the incident, transmitted 
and reflected wave packets move with the different velocities. As a rule, the transmitted 
packet moves faster than the incident one (for example, it takes place for an opaque 
rectangular barrier). In this case the CM of the transmitted packet appears behind of the 
barrier well before the incident packet arrives at the barrier region, i.e., before forming 
the reflected packet. The corresponding phase times were found to depend on the initial 
distance between the barrier and wave packet: the farther from the barrier is at the initial 
time the CM of the incident packet, the earlier the CM of the transmitted packet appears 
behind of the barrier. 

Thus, according to the SWPA, in the general case the transmitted particles seem to be 
accelerated (on the average) in tunneling, and the corresponding delay time (positive or 
negative) depends strongly on the initial distance between the particle and barrier. Note 
that the observed behavior of wave packets is entirely in accordance with the foundations 
of quantum mechanics, since it follows from the OSE to describe the tunneling process. At 
the same time, its interpretation mentioned above appears to be mistaken. It is intuitively 
evident that a static potential barrier must not accelerate (on the average) a particle, 
wherever it moves after the scattering event. Besides, the delay time must not depend on 
the initial distance between the particle and barrier. So that the behavior of wave packets 
in tunneling requires a reasonable explanation. The lack of such an explanation as well 
as the lack of well-defined characteristic times have a common root in the SWPA. 

At first sight, all the above means that the timing procedure for an everywhere free 
particle cannot be, in principle, used for a tunneling particle. However, in our opinion, 
this is not the case. One has only to bear in mind that the incident, transmitted and 
reflected wave packets describe quantum ensembles with the different number of particles. 
This fact have not been taken into account in the SWPA. The phase time for transmitted 
particles was deflned in this approach as the difference between the time of arrival of the 
CM of the transmitted packet at a point far behind of the barrier, and that of the CM of 
the incident packet at a point far before the barrier. Several authors have been cast doubt 
on the correctness of such a step. For example, considering the behavior of the incident 
and transmitted wave packets in the numerical modeling of tunneling, the authors of 
[14, 15] emphasize that "... arriving peaks do not turn into transmitted peaks In 
[1, 4, 6] it is pointed to "... the desirabihty of constructing an "effective" incident packet 
which would play the role of the counterpart to the transmitted packet..." Moreover, as 
was pointed out in [13], classical mechanics"... would suggest that the energy distribution 
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of the transmitted and reflected particles is the same before and after the coUision with 
the barrier..." 

To take into account these remarks, one needs to decompose the incident wave packet 
into (to-be-transmitted and to-be-reflected) parts. However, there have been an opinion 
that this step would contradict the basic principles of quantum theory. For example, 
Jaworski and Wardlaw [13] pointed out "... that quantum mechanically there is no sense 
in speaking about transmitted and reflected particles before they are detected as such..." 
As a result, the above idea remains unrealized. All the known attempts to describe 
individually transmitted and reflected particles have been made beyond the idea of the 
individual description of these subensembles at the initial stage of scattering. For example, 
Steinberg [10] offered to introduce the so-called conditional probabilities (however, being 
complex, these quantities lead to the serious problem of their interpretation) . There is also 
a radical viewpoint (see, for example, [16]) that an individual description of transmitted 
and reflected particles is inadmissible in quantum mechanics. 

However, as will be shown below, in the frame of the statistical interpretation (SI) 
of quantum mechanics (see, for example, [17]), the desirable decomposition does not 
contradict quantum theory. Besides, there is another aspect of the tunneling problem 
that also requires such a decomposition. We have in mind making use of the term "two- 
channel scattering" in solving the tunneling time problem. So. according to the definition 
given in [18, 19], this elastic one-particle one- dimensional scattering must be treated as 
an one-channel scattering. At the same time, according to Hauge and St0vneng [1], the 
term "scattering channel" in this problem should be associated with "... distinct final 
states (transmission and refiection ...)". By this definition, tunneling is a two-channel 
scattering. 

We agree entirely with this viewpoint. The mere fact that investigators attempt to 
define characteristic times individually for a transmitted and refiected particle means that 
they distinguish in this task two scattering channels. However, one has to recognize that 
such a viewpoint needs a proper theoretical basis. Namely, one should to develop a for- 
malism to describe individually the quantum subensembles of transmitted and reflected 
particles both at the initial and flnal stages of tunneling. Strictly speaking, it is illegiti- 
mately to say, without such a formalism, about two scattering channels in this case. In 
our opinion, just ignoring this fact leads to the above paradoxes. 

Note that there is a variety of scattering problems related to the electron transport 
(see, for example, [20]) in artiflcial quantum structures (e.g., in quantum wires) where a 
similar situation takes place. Right quantum wires which are nonuniform in some lim- 
ited spatial interval give the most simple generalization of the considered one-dimensional 
model. More complicated scattering problems arise for branching quantum wires where 
the number of scattering channels is larger than two. It is obvious that the above diffi- 
culties of the one-dimensional case must appear in these models too. To overcome them, 
one needs to treat the one-dimensional scattering as a two-channel one, and develop a 
relevant theory. 

So, in our opinion, the theory of tunneling must include an individual asymptotic 
description of the subensembles of transmitted and reflected particles. The purpose of 
our paper is to elaborate a desirable formalism in the framework of the SI (in our opinion. 
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solving the tunneling time problem is sensitive to the choice of some interpretation of 
quantum mechanics). Note that our ideas to solve this problem have been offered in 
papers [21, 22]. However some aspects of that solution must be corrected and developed 
further. In particular, we have to study here the role of spreading wave packets, which 
have been leaved out of account in [21, 22]. This effect is obvious to increase the scattering 
time. Moreover, when the velocity of the CM is smaller than that of the wave-packet's 
fronts, tunneling represents an incompleted scattering. 

The paper is organized as follows. In Section 1 we set the tunneling problem. In 
Section 2 we display explicitly the paradoxes arising in the standard approaches. For 
this purpose we analyze in detail the behavior of wave packets in the context of the 
SWPA. In Section 3 we show that the wave function describing the quantum ensemble 
of tunneling particles contains all information needed for the individual description of 
both scattering channels at the initial stage of scattering. We find here the corresponding 
counterparts for the transmitted and reflected packets, which describe individually both 
scattering channels at early times. In this section we answer also the question, how, in 
principle, to verify experimentally the formalism presented. In Section 4 we define delay 
times for transmitted and refiected particles. To estimate the duration of the scattering 
event, in Section 5 we determine the low bound of the scattering time. Besides, we 
derive here the condition which must be fulfilled for a completed scattering. In Section 
6, the properties of the characteristic times are illustrated in the case of the Gaussian 
wave packet tunneling through a rectangular barrier. Some useful expressions for the 
asymptotic expectation values of the position and wave-number operators as well as for 
their mean-square deviations are presented in Appendix. 

1 Setting the problem for a completed scattering 

Suppose that a particle moves from the left toward the time- independent potential barrier 
V{x) which confined to the finite spatial interval [a,b] {a > 0); d = b — a is the barrier 
width. Let the state of a particle be described at i = by the normalized wave function 
^^0(2^) belonging to the set S'oo which consists from infinitely differentiable functions van- 
ishing exponentially in the limit |a;| — 00. The Fourier-transforms of such functions are 
known to belong to the set 6*00 as well. This property guarantees that the position and mo- 
mentum operators both are well-defined. Let < \I'o|a;|\I'o >= 0, < ^'o|p|^'o >= ^ko > 0, 
Zq =< ^o|a;^|^E'o >; here Iq is the wave-packet's half-width a,t t — (Zq << a); x and p are 
the operators of the particle's position and momentum, respectively. 

Another important requirement should restrict the rate of spreading the incident wave 
packet. We must be sure that at early times all particles of the corresponding quantum 
ensemble move toward the barrier (only a negligible part of particles is assumed to move at 
t — away from the barrier). This means that the packet's spreading must be sufficiently 
ineffective (see condition (40) in Section 5). 

For our purposes it is useful also to introduce the following notations. Let H be the 
Hamiltonian of the problem: H = Hq + V{x), where Hq is the energy operator for a free 
particle. Besides, let Hq'^^ be the auxiliary Hamiltonian describing the refiection of a free 
particle off the absolutely opaque potential wall located at the middle point of the interval 
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[a,b], i.e., at Xmidp = (a + 6)/2. Let us introduce also the projection operators ©+ and 
©_ coinciding with the corresponding ^-functions: ©+ = 9{x — 6), ©_ = ^(a — x). 

In the general form, the solution of the corresponding temporal OSE can be written 
formally as e^^^*^^^ q{x) . To solve explicitly this equation, we will use the transfer matrix 
method (TMM) [23] that allows one to calculate the tunneling parameters for any system 
of potential barriers. Then, the (matched) solution describing the incident and reflected 
waves {x < a) with the given wave-number k can be written as 

^/e/t — [exp(iA;a;) + (prefi^) exp{—ikx)] exp[—iE{k)t/h], (1) 

where 



(j)ref{k) = \/R{k) exp 

the solution 



i{2ka + J{k)-F{k)-^) 



bright = K{k) exp[i(/cx - E{k)t/h)l (2) 
represents the transmitted wave {x > h)] 



{k) = ^jTiJ^)ex^[i{J{k) - kd)]] 



Here E{k) — ^^A;^/2m; T{k) (the real transmission coefficient) and J{k) (phase) are even 
and odd functions of /c, respectively; F(—k) = tt — F{k). 

Thus, for the temporal OSE, the solution to satisfy the above initial condition is given 

by 



"^leftix, t) = -j= I [finc{k, t) + fref{k, t)] exp{ikx)dk (3) 



(for the region x < a) and 

1 

'^right{x,t) = —7= ftr{k,t)exp{ikx)dk (4) 
V 27r J-oo 

(for X > b), where 

finc{k,t) = c ■ A{k; ko,lo) exp[-iE{k)t/h], 

fref{~k,t) — (f)ref{k)finc{k,t), 

ftr{k,t) = (t)tr{k)fincik,t); 

c is a normalization constant; the weight function A[k; /cq, Iq) is the Fourier-transform of 
^o(a;)- By the initial condition both ^o(a;) and A{k;kQ,lQ) belong to ^oo. In particular, 
for ^o{x) describing the Gaussian wave packet, A{k; ko, Iq) — exp {—ll{k — koY). 

Here functions (3) and (4) present wave packets moving in the out-of-barrier regions. 
Expression (4) describes the transmitted wave packet. Solution (3) consists of the in- 
cident and reflected packets (in fact, the latter appears only on arriving the incident 
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packet at the barrier region). To study these packets, it is convenient to pass into 
the /^-representation. It is obvious that in the case of a completed scattering, at late 
times, the functions ftr{k,t) and fref{k,t) approximate Fourier-transforms of the func- 
tions 0+e~*^*/^^o(3^) and Q^e~''^^^^'^Q{x), respectively. The function finc{k,t) approx- 
imates the Fourier-transform of e^''^^^^"^o(x) at early times. A detailed description of 
these packets in the k representation is presented in Appendix. 

As follows from Appendix (see (57), (59) and (60)), for a completed scattering the 
following normalization conditions must be valid (see also [13]). For early times 

/a roc 
\^ieft{x,t)\^dx ^ \finc{k,t)\^dk = l; (5) 
-oo J —oo 

for sufficiently large times 

/a poo _ 

\^left{x,t)\^dx^ / \fref{k,t)\^dk = R; (6) 

/ \^ri9ht{x,t)\^dx ^ / \ftrik,t)\^dk = f. (7) 
Jb J —oo 

Here T and R are the mean values of the transmission and reflection coefficients, respec- 
tively: f =< T{k) >inc', R =< R{k) >inc', angle brackets denote averaging over a relevant 
wave packet (see (55)). We will assume further that conditions (5)- (7) are satisfied. 

As was pointed out above, timing a particle can be made only beyond the scattering 
region. All needed information about the influence of the potential barrier on a particle 
is contained, via the tunneling parameters, in expressions (3) and (4). We will consider 
that these parameters have already been known. 

2 Paradoxes of the standard wave-packet analysis 

In order to display explicitly the basic shortcoming of the SWPA, let us rederive the 
tunneling times in the context of that approach. The only difference is that we use here 
the TMM [23]. 

As is shown in Appendix, the expectation values of the x-operator for all three packets 
are given by expressions (see (67) - (69)) 

< X >inc (t) = m'^hkot (8) 
for the sufficiently early times, and 

< X >tr (t) = m~^h < k >tr t + d- < J' >tr, (9) 
< X >ref (t) = -mr'^h < -k >ref t + 2a+ < J' - F' >ref (10) 

for the sufficiently large times. Note (see Appendix) that the prime denotes here the 
derivative with respect to k. 

Let Zi be a point to lie at some distance Li {Li ^ Iq and a — Li ^ Iq) from the left 
boundary of the barrier, and Z2 be a point to lie at some distance L2 {L2 ^ Iq) from 
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its right boundary. Thus, conditions (5) - (7) are fulfilled in this case. Note that the 
fulfilment of the condition a — Li ^ Iq must guarantee that (almost) all particles of the 
quantum ensemble are localized, at t = 0, to the left of the detector placed at the point 
Zi. In any case, a part of the quantum ensemble to violate this condition must be much 
smaller than T. 

Following the SWPA [1, 4, 5, 6], let us define the difference between the time of arrival 
of the CM of the incident packet at the point Zi, and that of the transmitted packet at 
the point Z2 (this time will be called below as the "transmission time"). Analogously, 
let the "reflection time" be the difference between the time of arrival of the CM of the 
incident packet at the point Zi and that of the reflected wave packet at the same point. 

Let ti and ^2 be such instants of time that 

< X >inc (tl) ^ a - Li, < X >tr {t2) ^ b + L2. (11) 

Considering (8) and (9), one can write then the "transmission time" Attr {^Ur — t2 — h) 
for the given interval in the form 



m 

Attr = - 

n 



< J' >tr +L2 Li / 1 1 

+ -; ha 



< k >tr ko \< k >tr ko 

For the reflected packet, let t[ and be such instances of time that 



(12) 



< X >inc (t'l) =< X >ref (4) = O " -^1- (13) 

From equations (8), (10) and (13) it follows that the "reflection time" Ai^e/ (Ai^e/ = 
^2 ~ ^'i) can be written as 



m 

Atref — 



< J' - F' >ref +L, Li / 1 



(14) 

< -k >ref ko \< -k >ref ko J \ 

Notice that the expectation values of k for all three wave packets coincide only in the 
limit Zo — ^ oo (i-e., for particles with a well-deflned momentum). In the general case these 
quantities are distinct (it is this asymptotic property of the wave packets that is usually 
treated [2] as the acceleration (or retardation) of a particle in tunneling). From Appendix 
(see (63)) it follows that the rule 

T- <k >tr +R- < -k >ref^ ko (15) 

must be true. 

As is seen, times (12) and (14) cannot serve as characteristic times for a particle. 
Due to the last term in (12) and (14), the above times depend essentially on the initial 
distance between the wave packet and barrier, with Li being fixed (a similar peculiarity 
arises in [13]). These contributions are dominant for the sufficiently large distance a. 
Moreover, one of these terms, either in (12) or in (14), must be negative. For example, 
it takes place for the transmitted wave packet in the case of the under-barrier tunneling 
through an opaque rectangular barrier. As was mentioned above, the numerical modelling 
of the tunneling process (see, for example, [1, 4, 5, 6]) shows, in this premature 
appearance of the CM of the transmitted packet behind of the barrier. This fact have 
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been treated in [2] as the evidence of the lack of a causal link between the transmitted and 
incident wave packets. To avoid this effect, the fields of application of the SWPA have 
been restricted (see [1, 4]) by the limiting case /q ~^ oo. A simple analysis shows however 
that the last terms in (12) and (14) remain dominant in the limit Iq oo, with the ratio 
lo/a being fixed. Note that the limit Iq ^ oo with a fixed value of a is unacceptable in this 
analysis, because it contradicts the initial condition a ^ lo for a completed scattering. 
Thus, even in the limit Iq ^ oo the above analysis carried out in the context of the SWPA 
does not provide the well-defined characteristic times for a particle. 

3 Formalism of the individual description of the sub- 
ensembles of transmitted and reflected particles at 
the initial stage of tunneling 

We think that a principle mistake made in the SWPA, in timing scattering particles, 
is that the "transmission time" was obtained there from the analysis of the relative motion 
of the transmitted wave packet with respect to the incident one. This is a physically mean- 
ingless step. The point is that the incident packet describes the whole quantum ensemble 
of particles, while the transmitted packet presents only its part. The incident packet can 
be used as a reference only for the transmitted and refiected packets taken jointly; in this 
case the initial and final states both describe the whole ensemble of particles. In the indi- 
vidual description of the transmitted (or reflected) packet, its motion should be compared 
with the corresponding counterpart which describes the state of transmitted (or reflected) 
particles at the initial stage of the scattering event. Searching for such counterparts for 
both subensembles is the following step of our analysis. 

Initial states of the subensembles of transmitted and reflected 
particles 

Note that /^-distributions for both subensembles at early times can be obtained imme- 
diately if one takes into account the physical sense of the transmission (and reflection) 
coefficient as well as the Born interpretation of a wave function. So, the transmission coef- 
ficient T{k) is, by definition, the probability for the incident particle with the momentum 
hk to pass ultimately through the barrier. The expression \finc{k, t)\'^dk is the probability 
for the particle to be, in k space, at time t in the interval [k,k + dk]. Since both the 
probabilities describe statistically independent events, the expression T{k)\finc{k,t)\'^dk 
is evident to give the probability that both these opportunities happen jointly for this 
particle (note that this joint probability does not represent a conditional probability in- 
troduced in [10]; the former is always positive unlike the latter). From this it follows that 
the function flnc{k.,t), where 

fl:,,{K t) = ^JT{k) ■ f,nc{k, t) eMiO{k)tr), (16) 
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may be considered as a wave function describing at the initial stage the subensemble of 
transmitted particles. Similarly, the function fl^{{k,t), where 



fZl{k,t) = ./R{k)-finc{k,t)exp{te{k)ref), (17) 

may be considered as a wave function to describe, at this stage, reflected particles. Here 
6tr{k) and 6ref{k) are arbitrary real functions. 

We see that the /^-distributions for the to-be-transmittcd and transmitted wave packets 
are the same. The similar situation takes place for reflected particles. The only difference 
is that the sign of k is different after the reflection. This property corresponds to that of a 
classical scattering pointed out in [13] (see Introduction). To complete the correspondence, 
one needs to state the initial condition for each scattering channel. Namely, we will 
consider that the CMs of the to-be-transmitted and to-be-reflected wave packets should 
start at t = from the same point, i.e, from the origin. 

One can see that the above conditions are not sufficient to find uniquely the phase 
shifts 9tr{k) and 9ref{k). It is easily to show that there is a variety of wave packets, with 
the different functions 9tr{k) and 9ref{k), which start from the origin and described by 
the same /c-distribution. However, for the following, of great importance is only the fact 
that the positions of CMs of all such packets coincide at any moment t. 

For example, the initial states of the to-be-transmitted and to-be-reflected subensem- 
bles, which obey the above requirements, can be deflned as follows 

*L(^) - ^L(^ - ^tr): ^:i{x) = ^Ziix - Xref), (18) 

where 

^Ux) = lim e'"°'/^e+e-'^'^^Mx), ^Indx) = lim e'^o''*/^e_e-'^*/^*o(a;); 

t—fOO t—^oo 

(19) 
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It is easy to show that the phases 9{k)tr and 9{k)ref of the Fourier transforms fi^dk, 0) 
and flncik, 0) of the initial states \l/-^g(a;) and "^l^iix) read as 

e{k)tr = J{k) -k< J'{k) >tr 

9{k)ref = J(k) - F(k) -k< J'{k) - F'{k) >ref ■ 

Functions (18) are evident to be non-orthogonal. In this case the incident wave packet 
can be presented in the form 

Mk, t) = fZcik, t) + fi:i{k, t) + fi:ii{k, t) (21) 

where 

= f 1 - sjT{k) ■ ^{iO{k)tr) 4W) ■ 0Xp{l9{k)ref) " ^k, t). 
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The non-orthogonality of the initial states of the transmitted and reflected subensem- 
bles of particles is connected ultimately to the fact that the future of a single particle 
from the incident quantum ensemble is unpredictable in quantum mechanics. The sorting 
of particles at the initial stage of scattering is, in fact, purely conditional. In other words, 
there is an interchange of particles between the subensembles corresponding to the to-be- 
transmitted and to-be-reflected wave packets. Nevertheless, in spite of this exchange, the 
number of particles in each subensemble remains constant. One can show that, in this 



case, the following relations are valid, 

^ find fine > + < finclfinc > = < finc\finc > • (22) 

< fZdfL > = < ftrlftr >, < fZi\fi:i > = < freflfref >; (23) 

< k >L=< k >tr; < k >Zi=< -k >ref; (24) 

< £ >Tnc it) = m-'h < k >L 4; (25) 

< X it) = m-'h < k >Zi -t. (26) 



As is seen, in (22) there are neither terms with nor terms to describe interfer- 
ence between the three contributions entering decomposition (21). Besides, we have to 
emphasize that, at this stage, there is no interference between the incident and reflected 
wave packets (the latter has not yet been formed). Thus, from these expressions it fol- 
lows that, long before the collision, the whole quantum ensemble of incident particles 
do indeed consists from two parts (see (22)). Relations (23) suggest that, at this stage, 
the to-be-transmitted and to-be-reflected parts are described by flnc{k,t) and flncik,t), 
respectively. For each subensemble, both the number of particles (see (23)) and the ab- 
solute value of the average wave- number k (see (24)) must be the same both before and 
after the collision. One can easily show also that 

f- < k >L +R- < k >\t= ko. (27) 

Besides, at t = the CMs of both the packets are located, in accordance with the initial 
conditions, at the point x — (see (25) and (26)). 

Thus, at the stage long before the scattering event, the incident wave packet can be 
divided into the to-be-transmitted and to-be-reflected ones whose evolution is described 
uniquely by relations (22) - (27). This provides solid grounds to say now that tunneling is 
a two-channel scattering, and (answering the questions in [2]), that the potential barrier 
does not accelerate (on the average) a particle. By our approach, the absolute value of the 
average momentum of particles must be conserved asymptotically not only for the whole 
quantum ensemble, but also for its transmitted and reflected parts. We have to stress 
that this conservation law was stated in the frame of conventional quantum mechanics, 
without new guesses. 



12 



On the experimental verification of the individual description of 
scattering channels in tunneling 

At first sight, it seems to be strange tliat, altliougli the need for an individual description 
of both scattering channels has been recognized (see Introduction) , a desirable formalism 
have not yet been elaborated. In our opinion, the main obstacle which arises here is the 
problem of the interpretation of quantum mechanics, rather than a purely mathematical 
one. For example, according to the orthodox interpretation, the above formalism contra- 
dicts the causality principle since the knowledge of the initial states for both scattering 
channels should mean that a particle "feels" the barrier yet before the interaction (see, 
for example, [13]). We disagree entirely with such a viewpoint. It itself contradicts Born's 
interpretation of a wave function. We think that a proper understanding of the above 
formalism can be attained only in the framework of the SI which is known to base en- 
tirely on the statistical. Born interpretation of a wave function. According to the SI, our 
formalism does not contradict the causality principle, and can be verified experimentally. 
Let us consider the last question in detail. 

In order to verify our formalism, it is sufficient to check the fact that the ensemble of 
incident particles can indeed be divided into two parts, with the expectation values of the 
particle's position changing in accordance with relations (25) and (26). To perform this 
check, one needs to obtain the momentum and position (further, in this section, these 
variables will be denoted by P and Q, respectively) distributions for the ensemble of 
incident particles. For the subensembles of transmitted and reflected particles, one needs 
to obtain only the P-distributions. All the above suggests carrying out the sufficiently 
large (strictly speaking, infinite) number of identical experiments to provide the needed P- 
and Q-data. In order to satisfy conditions (5) - (7) all measurements must be performed 
far enough from the barrier. 

The basic asymptotic property of the subensembles which can be used in this test- 
ing is that P-distributions for each scattering channel, before and after the scattering 
event, must be the same. Thus, one needs firstly to obtain experimental P-data for the 
subensembles of transmitted and reflected particles, to extract from them the correspond- 
ing P-distributions, and then to sort out incident particles into two subsets with the given 
P-distributions. The last step is the most important in this testing. Its peculiarity is that 
sorting the P-data should give two corresponding subsets of the Q-data. This means that 
the above P- and Q-data must represent the set of (P, (5)-pairs. It is naturally to suppose 
that such pairs can be obtained as a result of a simultaneous measuring of the particle's 
position and momentum. According to the SI, " ...there is no conflict with quantum theory 
in thinking of a particle as having deflnite [unpredictable] values of both position and mo- 
mentum..." [17]. Moreover, we have to add that the values of P and Q can be measured 
simultaneously (jointly) with an arbitrary accuracy (in our opinion, if one assumed that 
the Heisenberg uncertainty relation does not admit an accurate simultaneous measuring 
of these quantities, he would then to recognize that this relation cannot be verifled in 
quantum mechanics). We consider that, strictly speaking, there are no restrictions, in the 
conventional quantum mechanics, to measure jointly these quantities. In particular, this 
means that one can construct, in principle, such devices that will allow one to measure P 
and Q simultaneously, without the violation of the connection between the P- and Q-data 
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imposed by the Fourier-transformation. Thus, in hne with the terminology introduced 
in [24], in this question we adhere the second extreme view. In any case, the possibihty 
to measure jointly P and Q have been accepted, in principle, by many authors (see, for 
example, [17, 25, 26, 27]). 

So, let for the given instant t (P, Q)-pairs for incident particles, and P-data for scat- 
tered particles be available: it is supposed that the total number of the experimental 
data is sufficiently large and the (P, (5)-pairs have been accidentally numerated. So, one 
may consider the P-distributions for the transmitted and reflected subensembles have 
been known. Let also the P-scale be partitioned into sets of intervals with a sufficiently 
small width hAk. Then, let us consider for incident particles the (P, (5)-pairs with the 
P values to belong some interval [hk,h{k + AA;)], take from this interval the T'(A;)-th 
part of these pairs (for example, with the first numbers) and include this part into the 
to-be-transmitted subensemble. Note that the values of Q must not be taken into ac- 
count in sorting. Another part should be associated to the to-be-reflected subensemble. 
All the above should be done for all /i/c-intervals. Then, for incident particles of both 
scattering channels, we can calculate the expectation values of P and Q, for the given 
t, and check lastly expressions (22)-(26). Of course, in doing so, we wait that the P- 
and Q-distributions for each subensemble are connected by the Fourier-transformation. 
Otherwise, they do not correspond some wave functions. 

For the following it is important to consider another variant of testing. Namely, let 
all measurements be carried out at a given point Zi rather than at a given instant of 
time. Now we will measure the particle's momentum P and the value of parameter t at 
which the incident particle arrives at the given point Zi. Since before the interaction 
the incident particle is free, the P-shape of the corresponding incident wave packet must 
remain unaltered in time. This means that such measurements should give the same 
P-distribution for the incident packet, as in the above case. Thus, the only difference 
between both these variants is that now one needs to sort (P, i)-pairs instead of (P, Q)- 
pairs. After the mean values of t for each subensemble have been calculated expressions 
(22)-(26) can be checked. This variant is more suitable for testing characteristic times 
introduced below. 

4 Delay times for the subensembles of transmitted 
and reflected particles 

Now, for both scattering channels, we can offer new definitions of characteristic times 
to describe the influence of the barrier on a particle. For this purpose we have again to 
consider the particle's motion in a wide spatial interval containing the barrier region. In 
particular, let us calculate the (average) transmission time, Ttr, spent by a particle in the 
interval [^1,^2]. We have to take into account the fact that before the interaction the 
motion of the to-be-transmitted wave packet is described by the function e~*^°*/^^'*^^(a;) 
with the Fourier transform //^c(/c, t) (see expression (16), (18) - (20)). Then the time-of- 
arrivals ti and t2 to correspond the extreme points Zi and Z2, respectively, must obey 
the equations 
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< ^ >L (tT) ^a-L,; <x >tr (tf) ^b + L2. (28) 

Prom this it follows that 

Ttr{L,, L2) = tl' - if = -— (< J' >tr +Li + L2) . (29) 

Ti < K >tr 

Similarly, let the reflection time, r^^f , be the difference ^2"^-^ — t]^-^ where 

< X >Zi (tT^) =< X >ref (t^f) ^a-L,, (30) 
(remind that the to-be-reflected packet is described at the initial stage of scattering by 
e-'^ot/n^ref(^^-^ with the Fourier transform fZi{k,t) (see expression (17), (18) - (20)) 
One can easily show that 



(^1) = n<-k> ^< ^' ~ ^' ^^'^ ^^^'^ ■ ^^^^ 

re/ 

As was shown in [23, 28], the sign of the phase F' is opposite for waves impinging 
the barrier from the right. This case is similar to that when the wave moves from the 
left to the potential barrier Vinv{x) such that Vinv{x) = V{a + b — x). The corresponding 
reflection time rj:^^ can be written as 



m 



(+) 



m 



(^1) = (< J' + >ref +2Li) . (32) 

^ n < —K >ref 

One has to stress once more that these quantities cannot be treated, at L1 — L2 — 0, 
as the transmission and reflection times for the barrier region. There is no experiment to 
measure their values in this the disposition of devices at the boundaries of the 

barrier docs not provide a reliable identification of transmitted and reflected particles: 
conditions (5) - (7) are not fulfilled in this case. 

Note that times (30) - (32) are not suitable to describe the influence of the barrier on 
a particle, because they include the contributions of the out-of-barrier regions. However, 
they enables one to define time delays to characterize the relative motion of a scattered 
and corresponding free particles. To determine these quantities, we have to take into 
account the fact that beyond the scattering region the mean velocity of a free particle 
(taken as a reference for each scattering channel) should coincide with that of a scattered 
particle. Then for transmitted particles the delay time rf^i can be written, in terms of 
wave functions (19) and the tunneling parameters, as 



Similarly, for reflected particles the delay time r^g/ can be written as 

(_) _ m _ _ < > / < ^l^\p/m\^7nl > 



(34) 
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The delay time for the corresponding inverted barrier Vinv{x) is given by 




m 



(< J' + F' >ref -d) 



(35) 



h < -k > 



ref 



Expressions (33) and (34) (or (35) for the inverted barrier) should be considered as the 
definitions of delay times in our approach. As is seen, these times do not depend on the 
initial distance between the incident packet and barrier. Besides, they do not equal to 
zero for the 5-potential, since < J' >tr and < J' + F' >ref do not vanish in this case. 

It is evident that Xtr (see expressions (20)) can be treated as the spatial delay for 
the subensemble of transmitted particles, with Xtr =< J' >tr —d. Similarly, Xref with 

can be treated as the spatial delay for reflected particles. In the 
case of the inverted barrier, the spatial delay is equal to < J' + F' >ref —d. 

It is useful to compare respectively expressions (29) and (31) with (4.8) and (4.9) 
presented in [1]. Unlike our definitions, the integral in expression (4.8) is divergent at the 
point /c = if A(0; /cq; ^o) 7^ and T(0) ^ 0. For expression (4.9) a similar situation arises 
if A{Q- ko, lo) 7^ and R{0) ^ 0. See also [16, 29]. 

5 The low bound of the scattering time 

It should be noted that the above delays for both subensembles are accumulated during 
the interaction of a particle with the barrier. In this case the delay times themselves say 
nothing about the duration of this process. Thus, one has to define the third characteristic 
time of tunnehng, which will be further referred to as the low bound of the scattering 
time. It is obvious that this quantity cannot be defined individually for each scattering 
channel, because it should describe just the very stage of the one-dimensional scattering 
when a particle cannot be identified as incident, transmitted or reflected one. Besides, 
in this case one should keep in mind that the scattering event lasts until the probability 
to find a particle in the barrier region is noticeable. This means that in defining this 
quantity, one should to take into account spreading the corresponding wave packet. It is 
obvious that the scattering time docs not coincide with the time of staying the particle 
in the barrier region (see Introduction). Formally speaking, to define the scattering time, 
one needs to find such temporal interval for which conditions (5) - (7) are violated. 

So, let us define such instant of time, tstart, at which the distance between the CM of 
the incident packet and the left boundary of the barrier is equal to the half- width of this 
packet, i.e.. 



We will treat the instant t start as the time of starting the scattering event. We will assume 
that for t < tstart conditions (5) - (7) are fulfilled yet with a sufficient accuracy. 

As regards the end of the scattering event, one has to take into account that the 
transmitted and refiected particles move, on this stage, in the disjoint spatial regions, i.e., 
their wave packets do not interfere with each other. To define the corresponding instant of 
time tend, one has to use the quantities Str+ref and < (5,x)^ >tr+ref characterizing jointly 
both these packets (see (72) and (80) in Appendix). Here Str+ref is a total mean distance 



(a- < X >inc {tstart)f =< {(^x)"^ >inc {tstart)- 



(36) 
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between the CMs of the transmitted and reflected packets, and the corresponding nearest 
boundaries of the barrier; < [SxY >tr+ref is a total mean-square deviation of x averaged 

over the transmitted and reflected packets. 

Thus, let tend be the instant of time at which 

^tr+refitend) =< >tr+ref (tend)- (37) 

This deflnition suggests that for t > tend conditions (5) - (7) are fulflUed with a sufficient 
accuracy. 

Either of the above equations has two roots. A simple analysis shows that one should 
take the smallest root, in the case of equation (36), and the biggest root, in the case of 
(37). So, the searched- for solutions to (36) and (37) can be written in the form 

_m akp - ^/gfcg + (a2 - /g) < [Sky' >inc 

''''' ~ n ' ki- < (a-)2 ' ^"^^^ 

(remind that a ^ /q) ; 

_m bko-X + ^ll'^kl + - 2kohx + {P - /^) < {Skf >tr+ref 

" n ■ kl- < {5ky >,r+ref ^ ^ 

(see (81)-(84)). The low bound of the scattering time Tgcatt can be defined now as the 
difference tend - t start- 

A simple analysis shows that this quantity is strictly positive when the inequality 

{kl- < {5kf >,nc) > (40) 

is fulfilled. It should be considered as a condition for a completed scattering. It guaran- 
tees that (almost) all incident particles start a,t t — toward the barrier, and that the 
transmitted and reflected packets occupy, in the limit t ^ oo. the disjoint spatial regions. 
For a completed scattering the hierarchy tend > tgtart > is obvious to take place. The 
second inequality is guaranteed by (40). The first one must be valid, since the number of 
particles in the whole quantum ensemble is constant. Or, in other words, the quantum 
ensemble of particles cannot exit a region before entering it. In this case it is important 
to note that, in the hmit k^ — >< (Sk)^ >inc, 

m{a? — Iq) ma 

t start ~ 



2hkoa 2hko 

That is, expression (38) has no singularity in this limit. 

One can easily show that there is the optimal value of Iq at which Tscatt is minimal. 
The point is that, in the limit lo — > oo, the scattering time grows together with lo, and, 
at small values of this parameter, this time is large because of the fast spreading of the 
wave packet. If requirement (40) is violated, the transmitted and reficctcd packets must 
be overlapped at t ^ oo due to their spreading. As a result, the scattering event becomes 
incomplete. In this case the process of scattering a particle mimics that of escaping the 
particle from the barrier region. Our approach allows one to define the scattering time 
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for the incomplete scattering too. However this question is beyond the framework of our 
paper. 

In the limit Iq oo, expressions (38) and (39) are esensially simplified. In this case 
we have / ~ /q, /oq ^< {6k)^ >inc ■ Thus, the spreading effect may be neglected for narrow 
(in /c-space) wave packets. Taking account only of the dominant terms in (38) and (39), 
we obtain 

in 

Tscatt = ^ (2/o+ < J' >inc - < RF' >inc) ■ (41) 

For the corresponding inverted barrier the last term in (41) has an opposite sign. This 
term is nonzero only for asymmetrical potential barriers. 

One has to stress once again that the above definition yields the low bound of the 
scattering time. For it does not take into account the "tails" of wave packets. This means 
that, if T <^ 1 (or ^ <^ 1), the to-be-transmitted (or to-be-reflected) wave packet may 
occur in the scattering region when t ^ [t start ■i'tend\- To correct Tscatt for these cases, one 
needs to "reduce" the size of the packet's tails. Namely, for these cases a part of particles 
in the "tails" must be smaller than min(T', R). It is obvious that such a correction should 
increase the scattering time. However, we think this is superfluous. The point is that the 
quantity Tscatt defined above is sufficient to estimate the duration of the scattering process 
for the body of the quantum ensemble of particles. As regards its small part (transmitted 
or reflected), the instants of time and if" (or f-^^ and f^^) given by equations (28) (or 
(30)) provide a sufficient information about its motion. 



6 Tunneling the Gaussian wave packet through rec- 
tangular barriers 

To display some properties of tunneling, we have considered rectangular potential barriers, 
and investigated in detail the tunneling parameters of a particle whose initial state is 
described by the Gaussian wave packet (GWP). The weight function A{k-, ko, /q) in (3) is 
defined in this case by the expression 

A{k; ko, lo) = exp{-ll{k - /cq)^). 

As was pointed out above, in the general case the asymptotic value of the average 
wave-number of transmitted particles differs from that of reflected particles. One can 
show that for the GWP 



> 

<R' > 

< -k >ref= kQ + - 

4 

Supposing that 



< k >tr= ko + ... ^7 ; (42) 



< -k >ref= ko ■ (43) 



< k >tr= ko + {M)tr, < -k >ref = ^0 + {^k)ref, 
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we can rewrite relations (42) and (43) in the form 

f . {Ak)tr = -R ■ {Ak)ref = (44) 

Note that R' = —T'. The first equation in (44) coincides with rule (15) to be valid for 
any wave packet. 

Let us also derive several useful correlations for the mean-square deviations of the k 
operator. Using the relations (62) in Appendix, one can show that for the GWP 

< {Sk)'^ >inc— 



^.c,^2^ <T{k){k-kof >,nc (A^<r^>^ 3. 
< {5k) >tr^ — — + {Ak\,; 



^/c,^2^ _ < R{k){k-kof >inc {Ak)ref < R' >inc , ,^,^,2 
< [OK) >ref- ^ p ^ ^72 ^ p ^ ^ y^'^Jref- 

Further calculations (see (83)) yield 

< m' i, _ . (45) 

Now we dwell shortly on the numerical analysis of some peculiarities of tunneling (a 
more detailed numerical analysis is supposed to be done in the following). All calculations 
have been carried out on the basis of the TMM [23] (see also [22]). To clear up the 
role of the spatial localization of a tunneling particle, we have investigated the main 
features of the Zo-dependence of the tunneling parameters for the particular cases when 
V(x) = Vq = O.SeV, Eq — 0.02eV; m — O.OGTmg where mg is the electron mass. Fig.l 
shows the ratios of the mean wave numbers of the transmitted and reflected packets to 
that of the incident packet versus log{d/lo). As is seen, for particles with the well-defined 
momentum or position, the expectation values of the momentum for the transmitted 
and incident packets are the same. The reason is that in the case of the well-defined 
momentum the incident packet consists, in fact, of a single wave. In the second case the 
particle tunnels through the barrier without reflection. The point is that the contribution 
of high-energy harmonics into the GWP, for which the barrier is more transparent, grows 
together with d/lo. As a result, the transmission coefficient increases too (see Fig. 1). In 
the domain < \og{d/lo) < 2 a situation arises when the contributions of the transmitted 
and reflected waves are approximately equal. In this case it takes place the most distortion 
of the packet's shape as well as the maximal difference between the expectation values 
of the momentum for the transmitted and incident wave packet. In the case of a many- 
barrier structure, noticeable variations ofT{d/lo) can be observed if there are single-wave 
resonances near the point Eq {Eq — E{ko)). However, in any case, T — > 1, in the limit 
Iq — > 0, irrespective of the barrier's shape and value of ko- 
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Now we address to the spatial delays (J' — d)tr and {J' — d)ref- As was shown in 
[28], for the rectangular barrier of height Vq and width d, F' = and the derivative J' is 
determined by the expressions 

_ 2{k^ - k^)eKd + {e + k"^)^ smh{2Kd) 

«;[4A;2«;2 + {k^ + k^Y smh'^{Kd)] ' ^ ^ 



where k = ^^2m{Vo - E)/h, E < Vq; 

_ 2{k^ + k^)k^Kd - {e - K^f SlTl{2Kd) 

where k = ^j2m{E -Vo)/h^ , E > Vq. 

This enables one to explain the numerical data obtained for (J' — d)tr and (J' — d)ref- 
As is seen from Fig. 2, both quantities are equal only for particles with the well-defined 
momentum. It is important that in the limit Iq ^ the (average) spatial delay for a 
transmitted particle vanishes. This property, taking place for any barrier, is due to the 
fact that the average energy of a particle grows infinitely in this limit (such a particle 
passes freely through the barrier). As regards refiected particles, for rectangular barriers 
the spatial delay tends to [J'{k) — d]\k=o. 

It is interesting also to consider the case of the under-barrier tunneling, providing 
that the barrier's width grows but the wave-packet width is fixed. It is the very case 
which is usually analyzed in the literature (e.g., see [1, 13, 30, 31]) to demonstrate the 
" superluminal" propagation of particles in tunneling. Prom (46) it follows that for < Vq 
and Kd ^ 1, 

f ^ 2k-\ (48) 

i.e., in this case, for a particle with a well-defined momentum, this quantity does not 
depend on d. For E > Vq and nd^ 1, 

J' ^ + 

4k^K^ + ktpSin^Kd) ^ ' 

Besides, ior E — Vq 

^-3 4 + k^^d- ^ 



where ktop — yj2mVo/h^. From (48) - (50) it follows that for large enough d the function 
J'{k) is smooth ioi E < Vq and rapidly oscillates for > Vq. Besides, it grows sharply 
in going from the "under-barrier" domain to the "above-barrier" one. A simple analysis 
shows that the wider the barrier, the larger the contribution of waves with E > Voio the 
transmitted packet. In this case the transmission coefficient decreases until waves with 
E <Vq dominate in the transmitted wave packet. However, when these waves have been 
filtered off this packet, T and < k >tr does not depend on d (see Fig. 3 and Fig. 4). These 
quantities, as functions of d which was varied from 5 nm to 200 nm, have been calculated 
for Iq = 15 nm; other parameters are the same as for fig.l. As is seen, for the sufficiently 
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large d, < k >tr-~ hop- That is, the average energy of particles transmitted through a 
wide opaque rectangular barrier exceeds slightly the barrier height. Everywhere in the 
(i-interval investigated the average energy of reflected particles coincides approximately 
with that of incident particles; the point is that here R ^ 1. 

Besides, we found (see Fig. 5) that in the "under-barrier" domain < J' >jrfti< 
J' >ref^ 2.86 nm. In the "above-barrier" domain, the behavior of < J' >tr (d) is de- 
scribed qualitatively by expression (49). As regards reflected particles, < J' >ref amounts 
about 2.86 nm in both these domains. For the case considered, Tgcatt ~ 0.143 ps. 

Conclusion 

In this paper we have developed a new variant of the wave-packet analysis and proposed 
a solution to the tunneling time problem. Its key point is a formalism of the individual 
asymptotic description of transmitted and reflected particles at the stage preceding the 
scattering event. We have shown that this formalism does not contradict the principles 
of quantum mechanics, and it can be verified experimentally. According to our approach, 
the tunneling process, being an elastic one-dimensional scattering of one particle on the 
time-independent potential barrier, must be treated as a two-channel scattering. For each 
scattering channel, the momentum distribution must conserve asymptotically. 

It is shown that the above formalism gives a good grounds to solve the tunneling 
time problem. For each scattering channel we have defined delay times which describe 
the relative motion of the scattered and corresponding free particle moving, outside the 
scattering region, with the same average velocity. In addition, to characterize the duration 
of the interaction of a particle with the potential barrier, we define the low bound of the 
scattering time. It describes jointly both scattering channels. As is shown, there is an 
optimum initial width of the incident wave packet, leading to the minimal value of the 
scattering time. This time is strictly positive and describes the main part of particles of 
the quantum ensemble. 

All three characteristic times are defined for wave packets of an arbitrary width. They 
are obtained in terms of the expectation values of the position and momentum operators 
(in this sense, characteristic times in the given approach are quantities of a secondary 
importance). It is important that all integrals arising in our formalism have no singularity 
at A; = 0, however the transmission coefficient and function A{k; ko, Iq) depend on k in 
the vicinity oi k — 0. This property may be useful, in particular, in studying the "ultra 
Hartmann effect" (see [29]). 

As is shown, the scattering process is completed, if only the CM of the incident wave 
packet moves quickly enough (see condition (40). Otherwise, the problem of scattering a 
particle off the potential barrier turns, in fact, into that of escaping the particle from the 
barrier region. Note that it is quite possible to adapt our formalism to this case too. 
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Appendix: the asymptotic properties of a wave func- 
tion in the /^-representation 

The peculiarity of a completed scattering is that at the initial instant of time the incident 
packet is located entirely to the left of the barrier. After the scattering event there is 
two packets moving away from the barrier. In the limit t — > oo, they are located in the 
disjoint spatial regions. 

It is suitable to present the wave functions describing the incident, transmitted and 
reflected packets in the form 

*(a;,i) = ^/ dkf{k,t)e^'\ (51) 



/27r 

where f{k,t) e Soo] 

f{k,t) ^ M{k;ko,lo)eMiC{k,t)); 
M{k;ko,lo) and ^(/c, are the real functions. In particular, for the incident packet 

hk'^t 

Minc(k; ko, k) = cA(k; ko, k); ^inc(k, t) = — —. (52) 
For the transmitted and reflected packets we have 

Mtr(k; ko, k) = ^|T{k)Mir,c{k] ko, k); ^tr(k, t) = ^inc(k, t) + J(k) - kd; (53) 



TT 



Mref{-k, ko) = ^/R{k)Mi^^{k■, ko, lo)] Ufi'k, t) = ^inc{k, t) + 2ka + J{k) - F{k) 2 - 

(54) 

In the case of a completed scattering the above packets provide the asymptotic behavior 
of a wave function. 

For any Hermitian operator Q, Fourier transformation (51)-(54) enables one to deter- 
mine the evolution of the expectation value < Q >, 

at the stages preceding and following the scattering event, where ^ is one of the above 
wave packets. That is, both for the whole quantum ensemble and for its subensembles, 
we will calculate the expectation values of Q by means of the same rule (55). 

Strictly speaking, for the incident and reflected packets, the integrals in (55) should 
be calculated over the interval (— oo, a]. For the transmitted packet, one should integrate 
over the interval [6, oo). The point is that expressions (51) -(54) for these packets are 
valid only for the corresponding spatial region and corresponding stage of scattering. 
However, taking into account that the body of each packet is located, in the limit i — > oo 
or i — > — oo, in its "own" spatial region, we may extend the integration in (55) onto the 
whole OX-axis. Due to this step the description of these packets becomes very simple. 
At the same time, the mistake introduced in the formalism is expected to be negligible: 
the farther is the packet from the barrier at the initial time, the smaller is this mistake. 
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It vanishes in the hmits t — > oo. Thus, the asymptotic wave function for a completed 
scattering (for £ Soo) may be studied in the /^-representation. Now, on the basis of 
this representation, we can find the main characteristics of all three packets, which are 
desirable for the following. 

Normalization 

Note that 

1 poo poo 

< >= / dx dkdk'f*{k', t)f{k, t) exp[i{k - k')x] = 

ZTT J — oo J —oo 

/oo poo 
dk\f{k,t)\'^= / dkM\k;ko,lo). (56) 
-oo J—oo 

For each packet we have then the following norms. Since the particle is located at the 
initial time to the left of the barrier, we have 

/oo 
dkMl^{k-ko,k)^l. (57) 
-oo 

Then, allowing for (52), we have 

/oo 
dkA^ik^koJo). (58) 
-oo 

For the transmitted packet, 

/OO poo _ 

dkMlik; ko, k) = / dkT{k)Ml^{k; ko, lo) =< T{k) >,^,= f. (59) 
-oo J —oo 

For the reflected packet, 

/oo roo 
dkMlf{k; ko, lo) = / dkR(k)Ml^(-k; ko). 
-oo J —oo 

Having made an obvious change of variables, we obtain 

< /I/ >re/ = < R{k) >^nc= R- (60) 

From (57) - (60) it follows that 

f + R=l. 

The expectation values of the operators (n is the positive num- 
ber) 

Considering (51), one can find for all the packets that 

Q-qi roo 

< >= -i < *|— >= / dkM^(k;ko,lQ)k 

ox J-oo 
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(i.e., ^ is a multiplication operator in this case). Then for any value of n we have 

< *|A;"|* >=< /|A;"|/ >= / dkM^{k;ko,lo)k''. (61) 

J— oo 

Now we can treat the individual packets. From (61) and (53) it follows that 

< ftr\k"'\ftr > — < finc\T{k)k'^\finc > ■ 

In a similar way we find also that 

< fref\kVref >= (-1)" < finc\R{k)k^\finc >, 

and, hence, 

< T{k)k'' >inc^ f < A;" >tr, < R{k)k'' >inc^ (-1)"^ < k^^ >ref ■ (62) 

As a consequence, the next correlation is obvious to be valid 

< A;" >inc^ f<k^ >tr +R < {-kY >ref ■ (63) 

The expectation values of the operator x 

We begin again with the expressions to be common for all three packets. We have 

< >= — / dx dkdk'f*{k',t)f{k,t)xexp[i{k - k')x] (64) 

Substituting —i-^ exp{i{k — k')x) for the expression xexjp{i{k — k')x), and integrating in 
parts, we find that 

<%x\^>=i r dkr{Kt)^%^^ 

J —oo ok 

= , r ,kM(k;kM^^^^^ - r dUM^kMf-^. (65) 
J-OO dk J-oo ok 

Since the first term here is equal to 

'-M\k-MM\tZ = ^. 

we have 

< >^ - j^JkM\k-k,.kf-^ = - < ff-^\f > ■ (66) 

For the incident and transmitted packets, taking into account expressions (52) and 
(53) for ^{k,t), we obtain 

<i X ^inc — k ^inci (67) 

m 
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ht 

< X >tr= — < k >tr - < J'(k) >tr +d. (68) 

m 

Since the functions J'{k) and F'{k) are even, from (54) it follows that 

< X >ref= 2a+ < J\k) - F'{k) >ref -— < -k >ref ■ (69) 

Let, at the instant t, Str be the distance between the CM of the transmitted packet and the 

nearest boundary of the barrier, i.e., Str =< x >tr —b. Similarly, let Sref be the distance 
between the CM and the corresponding barrier's boundary for the reflected packet at the 
same instant: Sref = Oj~ < ^ >re}- From (68) and (69) it follows that 

Str{t) ^ — <k>tr-< J'{k) >tr -O, (70) 

m 

Srefit) ^ — <-k >ref - < J'{k) - F' (k) >ref -a. (71) 

Let us define now the average distance Str+ref{t) describing the both packets jointly: 

Str+ref{t) = T Str{t) + RSref{t)- (72) 

Considering (70), (71) and (63), we get 

Str+ref{t) = — < k >inc (73) 

m 

where b — a+ < J'{k) >inc — < R{k)F'{k) >inc (note that < J'{k) >inc— d, and 
< F'{k) >i„c= when V{x) = 0). 

The mean-square deviations in x-space 

Let us derive firstly the expression to be common for all packets. We have 

< >=7r dkdk'f*{k' ,t)f{k,t)x^ exp[i{k - k')x]. 

Substituting the expression — ^ exp{i{k—k')x) for that x'^ exp{i{k—k')x), and integrating 
in parts, we find 

d'f{k,t) 



/oo I 
dkr{k,t)- 
-oo 



dk^ 



Since 



d'f{k,t) 



dk'^ 

we have 



< 



/OO fOO 
dkM[M{i'f - M"] - i / dk[{M^y^' + M^C] (74) 
-oo J—oo 



(hereinafter the prime denotes the derivative with respect to k if the functions of two 
variables are written without the independent variables). One can easily show that the 
last integral in (74) is equal to zero. Therefore 
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/oo roo 
dkM\k;ko,lo)W{k,t)]^+ / dk[M'{k;ko,lo)f. (75) 
-oo J —oo 

Let, for any operator Q, < {6Q)^ > be the mean-square deviation < > — < Q >^; 
6Q = Q— < Q >. Then for the operator x we have 

< {Sxf >=< (In' Mf> + < {SC'f > . (76) 

Now we are ready to determine these quantities for each packet. Using (76) and 
expressions (52)-(54), one can show that for incident packet 

< {SxY >inc^< {WAY >inc +^ < {Skf >inc (77) 

(the first term here, in accordance with the initial condition, is equal to Zq); for the 
transmitted packet 

< {5xf >tr^< {l^Mtrf >tr + < {^J'f >tr " 

ht nh'^ 

-2- < {dJ'){dk) >tr +^ < {Skf >tr\ (78) 

for the reflected packet 

< {5xf >ref = < (In'Mreff >ref + < (Sf - 5F'f >ref + 

+2- < {5 J' - SF'){Sk) >ref +^ < {Skf >ref • (79) 

Let us determine now the mean-square value of (5x)^ averaged jointly over the trans- 
mitted and reflected packets: 

< {5xf >tr+ref^ f < {5xf >tr +R < (Sx)^ >ref (80) 

(note that < {Sx)^ >tr+ref¥'< {Sx)"^ >inc because < x >tr^< x >ref in the general case.) 
Considering (77)- (79), we can reduce this expression to the form 

ht hh'^ 

< {5xf >tr+ref^ - 2—% + ^ < i^kf >tr+ref, (81) 

where 

f = T < (In Mt,)' >tr +R < Qn'Mreff >ref 
+f < {5J'f >tr +R < {SJ' - 5F'f >ref, 

X = f < {6J'){6k) >tr +R < (5 J' - SF'){Sk) >,e/, (82) 

< (dky >tr+ref= f < {5kf >tr +R < {5kf >ref • (83) 
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The first two terms in the expression for P may be rewritten, using (53), (54) and the 
correlation T' + R' — 0, as 

f {In'Af >inc -\ < (ln'T)(ln'i?) >,„e + 

+f < {5j'f >tr +R < {6 J' - 5F'f >ref] (84) 

here the second term is positive and not singular both at the resonances and at the point 
k = 0. Since < (In'yl)^ >mc= ^l, we have > l^. 
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Figure captions 

Figure 1. The transmission coefficient (o) as well as the ratios < k >tr /ko (solid line) and 
< —k >ref /ko (dashed line) versus lg(d//o), for the rectangular barriers with Vq = 0.3 
eV, d = 5 nm; Eq = 0.02 eV, m = 0.067me; where rrie is the electron's mass. 

Figure 2. The spatial delays < J' — d >tr and < J' — d >ref foi' the transmitted (solid line) 
and reflected (dashed line) particles, respectively. The barrier's and particle's parameters 
are the same as for fig.l. 

Figure 3. ln(T) versus d, for Iq — 15 nm; the other parameters are the same as for fig.l. 
Figure 4. The ratios < k >tr /kg versus d, for the same parameters (see fig.3). 
Figure 5. The dependence of < J' >tr on d, for the same parameters (see fig.3). 
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